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ABSTRACT 

One of the central goals of multi- wavelength galaxy cluster cosmology is to unite all cluster observ- 
ables to form a consistent understanding of cluster mass. Here, we study the impact of systematic 
effects from optical cluster catalogs on stacked SZ signals. We show that the optically predicted 
F-decrement can vary by as much as 50% based on the current 2a systematic uncertainties in the 
observed mass-richness relationship. Mis-centering and impurities will suppress the SZ signal com- 
pared to expectations for a clean and perfectly centered optical sample, but to a lesser degree. We 
show that the level of these variations and suppression is dependent on the amount of systematics in 
the optical cluster catalogs. We also study X-ray luminosity-dependent sub-sampling of the optical 
catalog and find that it creates Malmquist bias increasing the observed F-decremcnt of the stacked 
signal. We show that the current Planck measurements of the F-decrement around SDSS optical 
clusters and their X-ray counterparts are consistent with expectations after accounting for the la 
optical systematic uncertainties using the Johnston mass richness relation. 

Subject headings: dark energy — galaxies: clusters: general — galaxies: clusters: intracluster medium 
— X-rays: galaxies: clusters — cosmic background radiation 



1. INTRODUCTION 

Measurements of the abundance of galaxy clusters as a 
function of their masses and redshift provides an impor- 
tant const raint on the nature of dark matter and dark en- 
ergy (e.g-lVikhlinin et ai]l20"o1l iVanderlinde et alJl2010t 
iSehgal et al.l 1201 1T L Joint analysis of multi-wavelength 
observations, including opti cal cluster catalogs and their 
Sunyaev-Ze l'dovic h (SZ- iSunvaev fc Zeldovichl Il972t 
iBirkinshawj 119991 iCarlstrom et al.l 120021 ) counterparts, 
will help realize the full cosmological potentia l of galaxy 
clust ers ([Cunha et all 120091 IRozo et all 120091: IWu et all 
[2011 . 

Optical galaxy cluster surveys have identified thou- 
sands of clusters down to a mass limit of ~ 
10 14 Mp) (e.g., the max BCG catalog of SDSS clus- 
ters, [KoesteFeral] [20073), and millimeter wave surveys 
have discovered hundreds of cl u sters using the SZ effect 
(e.g.. IVanderlinde et al.l 120101: iWilliamson et all 120111 : 
iMarriage et alll201lHAde et all2011al ). albeit to a higher 
mass limit due to instrumental noise. Using these cat- 
alogs, researchers apply mass-observable relations to re- 
late the true underlying halo mass to observed properties, 
like the galaxy member count in the optical (richness, 
Ng a i or A200) or the SZ decrement (y or F500). 

One can probe scaling relations down to masses be- 
low the detection limit by stacking the signal around 
known clusters. For instance, stacking has been used 
to the great benefit o f mass calibration in weak lens- 
ing a nd X-ray studies ([Sheldon et al.ll2007t iRvkoff et all 
120081 ). Similarly, the SZ/X-ray cluster scaling-laws 
and pressure-profiles were evaluated by Komatsu et al. 
(2011) and Melin et al. (2011), who stacked the SZ signal 
from WMAP data around known optical/X-ray clusters. 
These joint optical/X-ray/SZ analyses allow researchers 
to take advantage of the large volumes and mass ranges 
from optical cluster catalogs in combination with the 



lower sca tter in the mass observable relation in X-ray/SZ 
catalogs dShaw et al. 2008; Nagai ct al. 2006; Rasi a et all 
[20lllMotl et al II2005I) . 

The SZ signal recovered from stacking Planck data at 
positions of the maxBCG (jKoester et al.|[2007af) clusters 
shows a deficit of SZ signal compared to what is ex- 
pected from current m ass-richness scaling relationships 
(jAghanim et alll2011al ): t his discrepancy has been con- 
firmed using WMAP data ([Draper et alll2011f ). This dis- 
crepancy ma nifests itself differently for two mass-ric hness 
calibrations ([Johnston et al 2007; R ozo et all2009D both 
of which are based on the iSheldon et all ([2007D stacked 
weak-lensing mass measurements o f the maxBCG clus- 
ters. For the I Johnston et all (|2007l ) calibration, a simple 
reduction in the global weak-lensing mass ca libration by 
25% w ould eliminate the discrepancy. The IRozo et "all 
((20091) mass calibration requires a larger correction and 
a scalin g law that is not self-similar. Aghani m et all 
(j2011a| ) also show that a subset of the maxBCG clus- 
ters wit h measured X-ray lum inosities from the MCXC 
catalog (jPiffaretti et al.H2011f ) can match the predicted 
Y500 vs. richness scaling relationship, although they did 
not consider selection effects inherent in such a hybrid 
catalog; 

The lAghanim et all (|2011al ) analysis was based on a 
comparison of the observed Y500 around maxBCG clus- 
ters to two models with different mass-richness calibra- 
tions and without including optical systematics. They 
evaluated the impact of impurities in the optical cat- 
alog as well as scatter in the mass-richness relations 
and concluded that neither could account for the ob- 
served discrepancy ind ividually. Here, we broaden the 
lAghanim et all (j2011al ) analysis to include the uncertain- 
ties in the mass calibrations as well as the combined sys- 
tematic effects in optical cluster catalogs. Instead of two 
model predictions to compare against the data, we look 
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at the family of predictions which come from uncertain- 
ties in the calibrations and the ranges of systematics in 
optical cluster catalogs. 

There are numerous systematic effects in optical galaxy 
cluster catalogs. These include the cluster selection (as 
a function of mass M and rcdshift z) which comprises: 
completeness — the probability that a true halo will be 
detected; and purity — the probability that a detec- 
tion correct l y ide ntifies a halo rather than noise (e.g., 
IMiller et al.ll2005h . Cluster redshifts estimated using 
photometric data are uncertain which introduces scat- 
ter in the observed redshift. There is uncertainty in the 
mass-richness calibration as well as scatter. Finally, mis- 
identification of BCGs in the maxBCG cluster-finder pro- 
duces a ngular offsets between true and recovered cluster 
centers ([Johnston et al.ll2007| ) called mis-centering. Cen- 
tering o ffsets driven by other m echanisms (e.g., astro- 
physical: iSandersorreFal] [2Q03) are smaller than those 
caused by the BCG mis-identification, and so we do not 
consider them in this study. 

Using mock clusters taken from N-body simulations, 
we directly manipulate the purity, mass-scatter, scaling 
calibrations and their uncertainties. We then re-create 
the Planck richness stacking technique on these mock 
catalogs to create model Ysoo-richness relations and com- 
pare to the Planck observations. In fj^l we describe the N- 
body simulations and the suite of simulated optical clus- 
ter catalogs with various systematics, the mock Planck 
SZ observations, the mock X-ray observations and the 
stacking procedure. We then show the results of stack- 
ing the SZ signal for each systematic to explore how each 
systematic can individually affect the SZ signal (|j3J), and 
we compare to the Planck joint SZ-optical and X-ray 
analyses ( fl3.1j) . Throughout this paper, we assume a 
ACDM cosmology with a J1a=0.75 and i?o=0-71 unless 
otherwise noted. 

2. SIMULATIONS 

We begin with a simulated mass function and halo po- 
sitions from an N-body lightcone. We then impose ob- 
servables and realistic systematic effects to produce mock 
optical catalogs and then dress the halos with gas and 
simulate Planck SZ observations. 

2.1. N-body Lightcone 

To generate the mock SZ maps and galaxy catalogs, we 
begin with the halo positions from a large (N = 1260 3 
particles, 1000 [Mpc h^ 1 } 3 ) cosmological dark matter 
simulation. Cosmological parameters were chosen to 
be consistent with those measur ed from the five-year 
WMAP data (jDunklev et al.ll2009D combined with large- 
scale structure observations, namely eg = 0.8, Qm = 
0.264 and ilb = 0.044. The simula tion was carried 
out us ing the tree-particle-mesh code of lBode fc Ostrikerl 
(|2003l) . In total, the lightcone covers a single octant on 
the sky (~ 5000 deg 2 ) to a redshift of 3, containing halos 
with masses M FO f > 3 x 10 13 h^ 1 M . 

The simulation does not provide any observables (e.g., 
richness or SZ/X-ray luminosity). We do not use the 
halo masses output from this particle simulation directly 
but rather use the procedure described in Section 12.21 
With the mass resolution available from this simulation 
we can reproduce the properties of the maxBCG cat- 



alog, including systematics, for clusters with M500 > 
6 x 10 13 h' 1 M Q or iV 20 o > 20. 

2.2. Simulated Halo Catalogs 

In the Planck analysis (lAghanim et al.l l2011a| ). the 
clusters are binned according to their optical richness 
(-/V2007 Ncais)- Richness is defined as the number of 
bright red galaxies (within the E/SO ridgeline) insid e 
i?2oo that are brighter than 0.4 L* (|Koester et al.ll2007al) . 
Recall, that richness is an observed quantity and at any 
fixed value clusters can have a range of true masses (the 
mass scatter). 

The halo catalog provides a mass function and large- 
scale structure according to our chosen cosmology and 
similar to the observed universe. We cannot directly 
assign richnesses to these halo masses that match the 
observed scatter. Therefore we create a mock catalog 
of masses and richnesses and assign them to the N-body 
halos to preserve the large scale structure of the universe. 

For each richness we center a Gaussian probability dis- 
tribution function (PDF) in In(AIass). The center of the 
Gaussian is taken from a particular scaling relation; the 
width represents the scatter in In(Mass) at fixed rich- 
ness. We draw from these Gaussian PDFs to create a list 
of masses including scatter for each richness. We adjust 
the number of draws from each PDF to reproduce the 
halo mass function (e.g., we draw more times from the 
low mass bins). This provides a table of richnesses and 
associated masses with the same halo mass function as 
the N-body simulation. We sort the N-body halos and 
this table by mass. We associate the positions of the 
N-body halos to the drawn table based on this ordering. 
This produces our mock cluster catalog which includes 
large-scale structure and reproduces a particular choice 
of scaling relations and scatter. We use these masses to 
create the SZ profiles (Section 12. 4j) and X-ray luminosi- 
ties (Section [275]) . 



2.3. Optical Cluster Catalogs 

Using the procedure described above we create mock 
catalogs that are modified as follows to include system- 
atic effects: 

1. Mass-richness Calibration: We varied the rich- 
nesses of the halos according t o equation 26 
(and a ssociated uncertainties) from I Johnston et al.l 
(j2007l) : 



(M2OO 1^200) — -^200 1 20 



2(H) 



20 



(1) 



M 200 |2o = (8.8 ± 0A stat ± l.l sys ) x 10 13 h- 1 M Q 

a N = 1.28 ±0.04 



or equation 4 from lRozo et al.l (|2009() : 



(M 500 \N 2 oo) 



0.71 X 10 14 h- 1 Mr: 



exp(B M |jv) 



40 



(2) 



a M]N = 1.06±0.08 stot ±0.08 sys 
B M \y = 0.95 ± 0.07 sta t ± 0A0 sys 

We look at one and two a deviations from these 
mass calibrations. Masses are converted from M200 
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to M500 assumin g an NFW profile and mass con- 

and arc rcla- 



p rome z 

centrations from lDuffv et aIT(|2008D . 
tive to the critical density. 

2. Completeness: We vary the fraction of halos in 
bins of redshift and mass. 

3. Purity: We add into the halo catalogs an addi- 
tional number of false halos in bins of mass and 
redshift to create samples with different purities. 
We either vary the purity as a constant with mass 
and re dshift or matc h the pu blished maxBCG pu- 
rity of iKoester et all (|2007al) . 

4. Redshifts: We scatter the true halo redshifts by 
normal distributions with varying widths as large 
as <7z = 0.05. 



Center Offsets: We offset the center for a 
fraction of the cluste rs according to Eq. 10 in 
Uohnston et al.l (|2007l ). For the offset clusters, the 
actual amount o f the offset is described by Eq 8 in 
See also Figures 4 and 5 in 
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Johnston et al 



( 2007 ) 



2007). 



6. Mass Scatter We vary the width of the log-normal 
distribution of masses at fixed richness. 

Realizations of mock optical cluster catalogs are cre- 
ated to investigate the impact of individual system- 
atic effects. These in clude maxBCG-like systematics 
(IKoester et al.ll2007allb1 ; Uohnston et~alll2007t IRozo et al.l 
I2009T) and more general systematics that are constant in 
redshift and mass (or richness). We also create catalogs 
combining maxBCG systematics to compare to data. In 
our maxBCG-like mocks, the fraction of incorrectly cen- 
tered clusters ranges from 12% in the highest richness 
bins to 39% in the lowest richness bins with a mean off- 
set of 0.6Mpc corresponding to 3' for a cluster at the 
mean redshift of z = 0.2; the completeness and purity 
arc >90% above M 50 o > 1 x 10 14 h^ 1 M Q and have an 
estimated uncertainty of 2.5%; the mass scatter is 0.45 ± 
0.10, similar to Rozo (a ln{M )\N 200 = 0.45±g;^ (95% CL) 
at N 200 w40). 

2.4. Mock SZ Sky Maps and the Stacked Signal 

The halo SZ signals are generated using a the rmal pres- 
sure profile suggested bv (jArnaud et al.ll2010f) and used 
in the P lanck maxBCG stacking a nalysis (|Aghanim et al.l 
l2011al) . iBonamente et al.l (|2011f) compares the pressure 
profiles of 25 massive relaxed clusters observed in X-ray 
and with the Sunyaev-Zerdo vich Array (SZA) and find 
that they agree well with the lArnaud et al.l (|2010f) pro- 
file up to i?5oo- We project the profile along the line- 
of-sight to produce a compton-Y profile, scaled to the 
appropriate size for each halo redshift. Mock Planck 
observations were created in each frequency band using 
the appropri ate beam sizes, i nstrument noise and pri- 
mary CMB (|Ade et al.ll2011af) temperature anisotropy. 
We concluded that the 143 GHz channel reproduced the 
dominant features of the multi-frequency analysis, and so 
we restricted our analysis solely to this channel, which 
has a beam size of 7.18 arcminutc FWHM and a noise of 
0.9 fj, K~ degree. 



At the position of each optical cluster, we extracted 
the integrated thermal SZ signal Y 500 from each SZ sky 
map using a ma tched filter (e.g., iHerranz et al.l 120021: 
Melin et al.ll2006f) with an Arnaud profile ( Arnaud et al.1 
20101 ). the size of which is inferred from cither the John- 
ston or Rozo richness-mass s caling relations (same as 
used in lAghanim et al.ll2011a|) . We stacked these match 
filtered signals in richness bins; then the amplitude is cal- 
ibrated by comparing the spherical Y500 of the halos with 
the amplitude in the stacked SZ signal in the absence of 
systematics. We found t hat including an i ntrinsic scatter 
of 25% in y 500 - A/500 (|Shaw et al.ll2008f ) did not affect 
our results beyond increasing statistical uncertainties in 
individual catalog realizations and so we did not include 
this additional scatter in the following analysis. 

2.5. maxBCG-MCXC Subsample 

The Planck team studied a subset of the maxBCG 
catalog whose positions were matched to within ^3 
arcminutes of X-ray clusters from the MCXC catalog 
(jPiffaretti et al.l f201 If ) . Starting with the masses of our 
simulated halos we assign X-r ay luminositie s (Lx ) and 
scatter according to Table 1 in I Arnaud et al.l (|2010l ). We 
reproduce the scatters in the Lx scaling r elations (at 
fixed m ass and richness) that are observed in lRozo et al.l 
(|2009f) where (J\^l)\m ranges from 0.5 at low mass to 
0.45 at high mass and <J\^l)\n IS a constant 0.85 at all 
richnesses. We also vary the input Lx and scatter to 
asses the sensitivity to those parameters. We then se- 
lect subsets of the simulated halos which have the same 
redshift and Lx distribution as the MCXC subsample. 
This allows us to reproduce the MCXC subsample with- 
out needing to characterize the exact selection function 
which is undoubtedly complex as this catalog is drawn 
from heterogeneous X-ray data. We also ensure that the 
mis-centering for this mock MCXC-maxBCG catalog is 
truncated at 3'. The maxBCG-MCXC mock catalogs 
need not have the same scatter in the mass-richness re- 
lation as we imprinted into the full maxBCG mock sam- 
ples. This is because we imprint the observed scatter 
from Arnaud et al. (2010) directly onto the full catalog 
and then draw a sub-sample. For the MCXC/maxBCG 
mock subsamples, <7i n (L)|jv drops to 0.70 and oWjvniATaoo 
drops to 0.40. 

2.6. Correlations of Observables 

The large scatter in true mass at fixed richness (see 
Section 12. 3p induces a correlation between the observed 
X-ray luminosity and SZ signal. We note that this cor- 
relation is distinct from a secondary correlation in the 
scatter of observables. Figure [T] shows the cluster den- 
sity as a function of Lx and Y500 in various richness 
bins. Within each richness bin there is a strong corre- 
lation between the two observables. This will be crucial 
for understanding the joint maxBCG-MCXC subsample 
discussed in Section |3~2] 

Our simulation pipeline does not create correlated scat- 
ters in the observables at fixed mass. Such correlations 
are expected due to co mmon substructue within clusters 
and projection effects (|White et al.l l2010h however they 
arc likely seco ndary effects and ar e beyond the scope of 
this work (see lAngulo et all (|2012t) ). 
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Fig. 1. — Correlation between simulated X-ray luminosities (Lx) 
and SZ signals (Y500) in the richness bins used in this work. The 
contours are drawn where the number of clusters is 20% as large 
as the number at the mean value of Lx and Y500 for each richness 
bin (the center of each contour). The gray vertical line illustrates 
the approximate limit in X-ray luminosity reached by some of the 
surve ys used in the construction of the MCXC l lPiffaretti et al.l 
120111) where they overlap with maxBCG. 



3. RESULTS 

In Figure [5J we compare the stacked Y500 in our family 
of mock cluster catalogs to a "perfe ct" cluster catalo g 
that has been calibrated according to lRozo et al.l (|2009t ). 
The "perfect" catalog uses a single calibration and does 
not contain any of the systematics we discuss in Section 
2.31 This is identical to the model the Planck t eam used 
to compare to the data (jAghanim et al-lfeOllaj ). In each 
panel, the solid black line shows the average ratio (over 
multiple mock realizations) for models which apply the 
fiducial maxBCG values for calibration, mis-centering, 
purity, and mass scatter individually (as described at the 
end of Section |2~3")) . The gray bands show the range of 
models using the 1 and 2 a uncertainties on those param- 
eters. Dotted- lines show more general models (e.g., 70% 
purity independen t of mass). We also sho w the Planck 
data presented in lAghanim et al.l (|2011af ). We do not 
show rcdshift scatter and completeness since we found 
them to have negligible effects at maxBCG levels. 

Systematic uncertainties (2er) in the mass-richness cal- 
ibration result in up to 50% range in the model Y500 
measurements. This is because the Y500 values from our 
perfect catalog are calculated from a single mass cali- 
bration, while the model Y500S are calculated using the 
masses drawn from the calibration including 1 and 2 a 
uncertainties. 

Mis-centering suppresses (biases low) the model Y500S 
over the entire mass range, with the largest effect at low 
mass (~ 25% suppression). This can be understood from 
the convolution of the Planck beam (~ T full width at 
half of maximum) and the centering offsets which are on 
average ~ 3' at the median redshift of the optical sam- 
ple. The offsets are large compared to the Planck beam, 
which blurs out the SZ-signal after the convolution. The 
impact of this effect increases to ~ 25% at low mass, since 
the maxBCG mis-centering fraction is mass dependent. 
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Fig. 2. — A comparison of the stacked I500 in our family of mock 
cluster catalogs to a sin gle "perfect" cluste r catalog that has been 
calibrated according to Rozo ct al. (2009). The solid black lines 
show the model with maxBCG-like systematics included (individ- 
ually). The gray bands show the range of models after wc include 
the 1 and 2 <r uncertainties on the individual optical systemat- 
ics in addition to statistical uncertainties. Gray dotted lines show 
more general models, while the blue lines in the bottom panel are 
specific to the maxBCG/MCXC sub-sample. The error bars are 
the Planck data. Uncertainty in the mass calibration is the dom- 
inant effect on the model predictions, however impurity and mis- 
ccntcring both bias the model predictions towards lower values of 
Ysoo- On the other hand, X-ray luminosity selected sub-samples 
(e.g., the MCXC) show highly biased Y500 predicted values (com- 
pared to a perfect optical catalog). Sec Figures \3\ and [4] for the 
combined effects of these systematics. 



Impurities suppress (biases low) the amplitude of the 
model Y500S by in troducing pure noise int o the SZ maps. 
As also noted bv lAghanim et ahl (|2011a| ). high levels of 
impurity would be required to explain the discrepancy 
with the data. Just as important, the weak-lcnsing cal- 
ibration of the mass-richness relation would also be af- 
fected by large impurities which would lead to an en- 
hancement in the mass-richness relation. Since Y500 ~ 
Ms , high impurities could even cause the observed SZ 
signal to be enhanced compared to the systematics-free 
case (something neither we nor Planck detect). Accurate 
modeling of the impact of impurities on Y500 requires sim- 
ulating its effect on the weak-lcnsing calibration of the 
optical catalog. 

Th e stated uncertainty in mass scatter (jRozo et al.l 
I2009D does not have a significant impact on the SZ sig- 
nal recovered using a maxBCG-like catalog (see the gray 
band in the bottom panel of Figure [5]). However, the 
same can not be said for the MCXC-likc subsample (blue 
lines in the same panel and see Section I2.5[) . The X- 
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Fig. 3. — The Planck data (error bars) com pared to the single 
perfect model used in lAgh anim et al.ll20lTal) (blue line) and to 
the range of models (gray bands) after jointly combining all of the 
individual systematic effects seen in Figure [2] The naive perfect 
model predicts higher (on average) Y500 values compared to the 
models which include catalog systematics. The data are consis- 
tent with our model predictions within 1 a for the Johnston mass 
calibration. 

ray selection causes a Malmquist bias in low richness 
bins where the X-ray sub-sample preferentially contains 
brighter (and thus more massive) clusters. Figured] illus- 
trates that selecting clusters above some Lx limit (like 
the example shown by the gray line) preferentially selects 
clusters with high 5*500 • Larger mass scatter increases the 
correlation between Lx and Y500 and therefore enhances 
the Malmquist bias. Richness bins that lie completely to 
the right of the Lx limit arc not affected by this bias and 
so the SZ signal there is not enhanced. 

3.1. Simulating Planck-maxBCG Joint Analysis 

Figure [3] compares the Planck results to our models. 
The Pla nck data (error bars) are the same in both pan- 
els from lAghanim et al.l (|2011aD . The solid blue lines 
shows the single naive perfect model based on either the 
Johnston (left) or Roz o (right) mass calibratio n in the ab- 
sence of systematics (jAghanim et aLlfeOllal ). The gray 
bands show model predictions based on our Monte-Carlo 
mock cluster catalog realizations which include all of the 
maxBCG optical catalog properties, uncertainties, and 
systematics shown in Figure [2] and which were applied 
in the original weak-lensing richness mass calibrations. 
While the Planck data are statistically inconsistent with 
the naive perfect model prediction, they lie at the lower 
edge of the models which include the ~1 a systematic 
uncertainties for the Johnston mass calibration. 

3.2. Simulating maxBCG-MCXC Joint Sample 

Figure [4] shows our prediction for the MCXC sub- 
sample of the maxBCG catalog compared to the Planck 
data. The gray bands here include simulated optical and 
X-Ray systematics as well as the X-ray selection func- 
tion. As expected from Figure [5}bottom, we see a bias 
in the predicted Y500 with decreasing richness due to the 
Malmquist bias present in low richness bins after only 
the brightest Lx are selected (see Figure [lj . The Planck 
observations lie inside the lower edge of the models which 
include the 1 and 2er systematic uncertainties for John- 
ston and Rozo mass calibrations respectively. 



4. DISCUSSION 
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Fig. 4.— The Planck data for the maxBCG/MCXC X-ray sub- 
sample (error bars) co mpared to the single perfect model used in 
(Aghanim et al. 2011a) (blue line) and to the range of models (gray 
bands) after jointly combining all of the individual systematic ef- 
fects seen in Figure [2] While the perfect model is the same as 
in Figure \3\ the gray bands here include the bias seen in Figure 
[2] (bottom) , which is caused after sub-sampling clusters based on 
their X-ray luminosities to match the observed data. The naive 
perfect model predicts lower (on average) Y500 values compared to 
the models which include catalog systematics. The data are consis- 
tent with our model predictions at the 1 (2) <r levels on the optical 
systematics for the Johnston (Rozo) mass calibration. 



The Planck team reported that the stacked SZ signal 
around optical clusters lies well below the single model 
expectation which does not include the optical catalog 
systematic uncertainties. On the other-hand, they find 
that the observed stacked Y500 values around an X-ray 
limited sub-sample are consistent with the naive optical 
model. They concluded that the gas properties of clus- 
ters appear to be more stably related to each other than 
the ga s-to-optical properties of clusters (jAghanim et al.l 
l2011a[ ). In this work, we reach a fundamentally differ- 
ent conclusion: the Y500 values observed by Planck are 
consistent with the model predictions for both the entire 
cluster sample and the X-ray sub-sample to within the 1 
a opti cal systematic uncertainties of the Uohnston et al.l 
mass calibration. Not only do we argue that there 
is no significant discrepancy between the models and 
the observed Planck stacked Y500 values around optical 
clusters, but we also argue that the optical and X-ray 
selected sub-samples simultaneously agree with model 
predictions. For instance, we can apply a single mass- 
richness calibration to the data and fit the predicted Y500 
models in Figures [3] and [4] simultaneously. However we 
do not pursue a joint SZ-optical mass calibration here, 
as it beyond the scope of this work. 

We find that the dominant source of optical systematic 
uncertainty comes from the mass calibration, which alone 
can account for most of the original discrepancy noted by 
lAghanim et alJ (|201 lal) . Impurities and centering errors 
combine to bias the model predictions towards lower Y500 
for the optical samples while mass scatter biases the pre- 
dictions high for low richness systems in the X-ray limited 
subsample. When fully accounted for, these systematics 
allow for models which are matched by the observed data 
for both the optical and X-ray cluster sub-samples in the 
Planck data. The range on the acceptable models is quite 
large and we note that the SZ-optical scaling laws can- 
not by precisely characterized using this type of stacking 
until the optical systematics improve (specifically mass 
calibration and its scatter). 

This work highlights the importance of multi- 
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wavelength studies of cluster properties as a source of 
cross-checks and a calibration. It is clear that optical sys- 
tematics cannot be ignored and future analysis of stacked 
clusters should be done using Monte Carlo analysis to in- 
clude a larger suite of systematic errors. 
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